UNCLASSIFIED 


AD  JOS  197 


DEFENSE  DOCUMENTATION  CENTER 

FOR 

SCIENTIFIC  m  fECHNICAL  INFORMATION 

CAMERON  STATION.  ALEXANDRIA.  VIRGINIA 


UNCLASSIFIED 


NOTICE:  When  goveraoent  or  other  drawings,  sped* 
flcatlons  or  other  data  are  used  for  any  puzpose 
other  than  In  connection  with  a  definitely  related 
goveznoient  procurement  operation,  the  U.  S. 
Government  thereby  Incurs  no  responsibility,  nor  any 
obligation  i^tsoeverj  and  the  fact  that  the  Govezn- 
inent  may  have  fozsulated,  fuznlshed,  or  In  any  way 
supplied  the  scdd  drawings,  specifications,  or  other 
data  Is  not  to  be  regarded  by  Implication  or  other¬ 
wise  as  In  any  manner  licensing  the  holder  or  any 
other  person  or  cozporatlon,  or  conveying  any  rl^ts 
or  pemlsslon  to  manufacture,  use  or  sell  any 
patented  Invention  that  may  In  any  way  be  related 
thereto. 


406  197 


(5)3- J-  6 


R  945 


DAVIDSON 

LABORATORY 


Report  9^5 


METHOD  OF  PREDICTING 
COURSE  STABILITY  AND 
TURNING  QUALITIES  OP  SHIPS 


by 

Wlnnlfred  R.  Jacobs 


R  9^5 


DAVIDSON  LABORATORY 
REPORT  9^5 


March  1963 


METHOD  OF  PREDICTING  COURSE  STABILITY 
AND  TURNING  QUALITIES  OF  SHIPS 


by 

Wlnnlfred  R.  Jacobs 


Prepared  under  sponsorship  of 
Office  of  Naval  Research 
Contract  Nonr  263(  ^O) 
Technically  administered  by- 
David  Taylor  Model  Basin 
DL  Project  LZ2388/04l 


Reproduction  In  whole  or  in  part  Is  permitted 
for  any  purpose  of  the  United  States  Government 


Approved 


I 


ABSTRACT 


An  analytical  method  combining  simplified  potential 
flow  theory  and  low  aspect-ratio  wing  theory  with  empirical 
modifications  for  a  real  viscous  fluid  Is  used  to  predict  the 
stability  derivatives  (first  order  hydrodynamic  force  and 
moment  derivatives)  of  a  family  of  hulls  In  order  to  estimate 
the  dependence  on  geometric  characteristics  of  course  stability 
and  turning  or  steering  qualities*  The  hulls  are  Taylor 
Standard  Series  forms  with  after  deadwood  removed  and  have 
the  same  length  and  prismatic  coefficient  but  varying  length- 
draft  and  beaim-draft  ratios  and  skeg  area.  Comparison  between 
the  values  calculated  by  this  method  and  those  obtained  from 
experimental  measurements  shows  good  agreement.  The  analytical 
method  can  predict  the  relative  effects  of  the  geometrical 
characteristics.  Calculated  magnitudes  are  slightly  different 
from  the  experimental  but  are  on  the  conservative  side.  How¬ 
ever,  since  the  hulls  tested  have  the  same  prismatic,  the  em¬ 
pirical  modification  for  the  rotary  moment  derivative  which  Is 
a  function  of  prismatic  coefficient  has  not  been  fully  tested. 
Necessary  refinements  of  the  method  must  wait  on  analysis  of 
data  on  series  of  forms  with  different  prismatic  as  well  as 
length-draft  and  beam-draft  ratios. 
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NOMENCLATURE 


A 

AR 

B 


D 


^6 


F 


g 

H 

h 

hfp 


|u^tH 


ki 

kg 

k» 

L 


profile  area  of  wing  or  hull,  ft^ 
aspect  ratio  of  wing 
beam,  ft 

dimensionless  lift  coefficient  based  on 
profile  area 

two-dimensional  lateral  added  mass  coef¬ 
ficient  (sectional  inertia  coefficient) 

diameter  of  turning  circle,  ft 

total  resistance  coefficient  of  the  hull 

force,  lb 

measured  lateral  force  coefficient 

acceleration  of  gravity 
maximum  draft,  ft 
local  draft,  ft 
maximum  skeg  height,  ft 

moment  of  Inertia  of  hull,  Ib-ft-sec^ 

added  moment  of  inertia  of  entrained  water 
(see  text),  Ib-ft-sec^ 


Lamb ^3  coefficients  of  accession  to  Inertia, 
longitudinal,  lateral  and  rotational 

lift,  lb 

dimensionless  lift  coefficient  based  on 
area  ^  •  H 

length,  ft 

mass  of  hull,  slugs 
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hull  mass  coefficient 


ml  = 


m[  = 


m 


mjj,  =  +  m» 

=  m^  +  m^ 


"•z 

N 

N»  = 

n’  = 


N 


^0+ 


H 


R 

Rf 

«r 

r»  = 


o  -  Ut 

3  -  X 


t 

u 


u  =  U  cos  /3 
u  =  -U  sin  6 

y 

X,  y,  Z 


longitudinal  added  mass  coefficient 

lateral  added  mass  coefficient  ( see  text) 
longitudinal  virtual  mass  coefficient 

lateral  virtual  mass  coefficient 
rotational  added  mass  coefficient  ( see  text) 
yawing  moment,  Ib-ft 

dimensionless  yawing  moment  coefficient 


virtual  moment  of  inertia  coefficient 

radius  of  turning  circle,  ft 
frictional  resistance,  lb 

residual  resistance,  lb 

dimensionless  angular  velocity 

dimensionless  distance  along  the  path  of 
the  center  of  gravity  of  the  hull 

time,  seconds 

velocity  of  the  center  of  gravity  of  the 
hull,  ft/sec 

x-component  of  U,  ft/sec 

y- component  of  U,  ft/sec 

coordinate  axes  fixed  in  the  hull  with 
origin  at  the  center  of  gravity 

distance  from  LOG  of  center  of  gravity  of 
lateral  added  mass,  ft 

distance  from  LOG  of  center  of  pressure  of 
lateral  force  Y,  ft 

distance  from  LOG  of  center  of  pressure  of 
tall  surface  or  skeg,  ft 
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Xg,  x-coordlnates  of  stern  and  bow,  respectively 

y  lateral  hydrodynamic  force,  lb 


i3 

6 

A 

P 


lateral  hydrodynamic  force  coefficient 

yaw  angle  or  drift  angle 
rudder  angle 

displacement  of  hull,  lb 

mass  density  of  the  fluid,  slugs/ft^ 

stability  indices 


Subscripts  (other  than  those  in  above  definitions) 


H 

refers 

to 

bare  hull 

1 

refers 

to 

ideal  fluid 

r» 

refers 

to 

derivative  with 

respect 

to  r* 

T 

refers 

to 

tail  surface  or 

skeg 

refers 

to 

derivative  with 

respect 

to  jS 
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INTRODUCTION 


The  object  of  the  program  whose  results  are  re¬ 
ported  here  was  to  develop  an  analytic  method  for  estimating 
the  course  stability  and  turning  or  steering  qualities  of 
ships.  The  method,  following  Martin,^  was  to  be  based  on  a 
combination  of  simplified  flow  theory  and  low  aspect-ratio 
wing  theory  and  was  to  be  assayed  by  comparison  with  signifi¬ 
cant  empirical  results. 

Measurements  exist  of  lateral  force  and  moment,  on 
straight  course  and  in  turn  (rotating  arm  tests),  on  a  series 
of  eight  models  (the  "84o"  Series)  having  the  same  parent  as 
the  Taylor  Standard  Series  without  the  deadwood  (faired- in 
skeg)  aft  (Fig.  l) .  The  models  were  of  the  same  length  and 
prismatic  but  with  varying  beam  and/or  draft  and  hence  dis¬ 
placement.  In  addition,  there  are  measurements  taken  on  the 
models  with  various  flat  plate  skegs  added  (Fig.  2),  but  no 
other  appendages. 

Data  measured  in  1946-47  were  reported  in  refs.  2 
and  3  which  were  concerned  primarily  with  equilibrium  turn¬ 
ing  conditions.  Measurements  were  limited  in  some  cases  to 
a  very  narrow  range  of  yaw  angle  around  the  equilibrium  turn¬ 
ing  angle  for  a  given  diameter  of  turn.  For  this  reason,  the 
data  on  models  with  skegs  reported  in  ref.  3  are  useless  for 
the  purpose  of  this  report.  However,  unpublished  experimental 
data  over  a  wide  range  of  yaw  angle  at  several  turning  radii, 
obtained  in  1951  at  this  laboratory,  are  available  for  two  of 
the  models  with  and  without  a  plate  skeg  extending  to  the 
stern.  A  third  hull  was  tested  here  in  1959^  without  skeg 
and  with  skegs  of  various  sizes,  and  an  analysis  of  the  data 
was  reported  by  Tsakonas.^ 

Martin^  recently  reanalyzed  the  straight  course 
data  obtained  in  1946-47  on  the  eight  bare  hulls  along  the 
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lines  suggested  by  Thieme,^  Inoue®  and  Fedyaevsky  and 
Sobolev.’^  His  method  of  treating  the  problem  of  ship  motions 
in  the  horizontal  plane  is  based  on  a  modified  low  aspect- 
ratio  wing  theory,  the  Munk  ideal  moment  and  cross-flow  drag 
theory.  The  last  was  used  in  estimating  the  nonlinear  force 
and  moment  range. 

The  present  work  considers  only  stability  deriva¬ 
tives  of  the  first  order  so  that  the  nonlinear  variations, 
the  quadratic  and  higher  terms,  may  be  neglected.  Also, 
since  in  the  linear  range  the  straight-course  data  are  much 
sparser  than  the  rotatlng-arm  data,  it  was  decided  to  give 
more  weight  to  the  latter  in  estimating  static  stability  de¬ 
rivatives.  It  has  been  the  experience  at  Davidson  Laboratory 
that  entirely  reliable  static  force  and  moment  rates  for 
straight- course  motion  can  be  obtained  from  rotatlng-arm  data 
at  sufficiently  large  turning  radii.  This  was  confirmed  in 
ref.  4  where  rotating-arm  results  were  compared  with  straight- 
course  data  measured  in  Tank  No.  3  of  Davidson  Laboratory 
using  the  same  towing  system  and  measuring  devices  as  for  the 
rotatlng-arm  experiments  in  Tank  No.  2.  In  this  way,  errors 
attributable  to  Inconsistent  mechanisms  were  avoided.  The 
earlier  straight-course  data  had  been  measured  in  Tank  No.  1 
using  different  towing  and  measuring  apparatus. 

Landweber  and  Johnson®  have  shown  that  the  simpler 
methods  derived  from  an  alliance  of  potential  flow  theory 
and  low  aspect-ratio  wing  theory  can  estimate  the  stability 
derivatives  as  accurately  as  sophisticated  methods  based  on 
more  realistic  theoretical  considerations.  The  more  complex 
methods  employed  in  ref.  8  require  as  many  simplifying 
assumptions  and  empirical  correction  factors  if  not  of  the 
same  kind. 


The  analytical  method  adopted  here  combines  the 
force  and  moment  approximations  which  were  derived  for 
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spheroids  from  potential  flow  theory  by  Lamb®  and  for  long 
slender  bodies  with  tapered  or  pointed  ends  by  Breslln^® 
with  Albrlng’s^^  empirical  modifications  for  viscid  flow. 
This  approximate  method  Is  based  on  simple  concepts,  yet 
correlation  can  be  considered  good  between  calculated  and 
experimentally  obtained  values.  Necessary  refinements  must 
wait  on  the  availability  of  further  data  on  hulls  of  other 
prismatic  with  and  without  skegs. 

This  project  was  sponsored  by  the  Office  of  Naval 
Research  under  Contract  Nonr  263(4o)  and  technically  ad¬ 
ministered  by  David  Taylor  Model  Basin. 


THE  ANALYTICAL  METHOD 


Simplified  Flow  Theory 


In  the  potential  flow  theory  the  hydrodynamic  force 
and  moment  rate  coefficients,  or  stability  derivatives,  of  an 
elongated  body  of  revolution  without  appendages  are  defined 
for  the  linearized  region  of  small  angles  of  attack  and  large 


radii  of  rotation  as:  on  straight  course. 


r'  = 


=  0 


=  =0 


■"H 


"H 


N' 


=  mi  -  m,’  =  N; 


(Munk  Ideal  moment) 


In  turn,  around  |S  =  0, 


■H 


=  0 


(l) 


•H 


The  notation  Is  that  of  the  Society  of  Naval  Architects  and 
Marine  Engineers  (see  Nomenclature  and  Pig.  3).  The  measured 
lateral  force  coefficient  Is  defined  as 

=  Y'  -  (m^  +  m’)  r' 
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and  its  derivative  with  respect  to  r*  as 


SF4 

-  (mi  +  m;) 


(2) 


where  is  the  mass  coefficient  of  the  hull  and  m^f  is  the 
longitudinal  added  mass  coefficient.  Lamb,®  considering  the 
added  mass  term  as  a  hydrodynamic  force,  defines 
y*  f  =  -  m{  =  -  k;,m’  where  k,  is  the  coefficient  of  longi- 

Xji  j-  ^  o 

tudinal  accession  to  inertia.  Equations  1  are  equivalent  to 
those  derived  by  Breslin^°  for  a  long  slender  body  with 
tapered  or  pointed  ends  from  three-dimensional  singularity 
distributions . 


The  first  of  eqs.  1  is  known  to  be  in  serious  error 
in  a  viscous  fluid.  To  quote  Arnsteln  and  Klemperer:  "When 
an  airship  is  propelled  at  an  angle  of  attack,  lift  forces 
are  created  in  a  similar  manner  as  by  the  wing  of  an  airplane. 
It  is  true  that  the  airship  *s  shape  as  a  wing  is  very  poor 
and  its  aspect  ratio  extremely  small;  but  the  size  of  the  ex¬ 
posed  surface  is  so  great  that  tremendous  aerodynamic  force 
components  at  right  angles  to  the  flight  path  can  be  evoked. 

Pedyaevsky  and  Sobolev^  have  defined  the  forces  and 
moments  acting  on  a  ship  by  identifying  the  body  of  the  ship 
with  a  wing.  In  this  analogy  the  span  of  the  wing  is  assumed 
to  be  double  the  draft  of  the  ship  to  take  into  account  the 
action  of  the  free  water  surface.  Tsakonas"^  shows,  by  a  com¬ 
parison  with  the  wind-tunnel  results  of  Flax  and  Lawrence, 
that  this  "solid  wall”  method  of  accounting  for  the  free  sur¬ 
face  effect  is  correct  for  moderate  speeds  when  the  Influence 
of  wave  making  may  be  neglected. 

Albring^^  has  derived  approximate  formulas  for  the 
stability  derivatives  of  a  body  of  revolution  with  and  without 
appendages  moving  in  a  viscous  and  eddying  fluid.  The  lift  L 
on  a  bare  hull  is  the  force  developed  by  the  "imagined  angle  of 
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of  attack  at  the  stern"  of  a  "correspondingly  shaped  solid 
without  the  effect  of  curvature."  It  is  assumed  or  taken 
from  experimental  measurements.  It  acts  at  a  distance  x 
from  the  center  of  gravity  and  Albrlng  suggests  that  for  most 
bodies  of  revolution  Xp  =  0  Is  a  good  approximation.  For  a 
bare  hull  his  results  are  given  In  ref.  11  as: 
on  straight  course,  r'  =  ^  =  0,  ^ 


L'  assumed  or  measured 


^6 


y  (3) 


NJ,  =  m2-mJ+-^L'  «mi-m/«NJ, 


In  turn,  around  |3  =  0 


A 


»  0 


■’H 


Albrlng  does  not  give  a  formula  for  N' ,  for  the 

body  without  appendages.  If  one  assumes  that  the  center  of 
pressure  does  not  shift  when  going  from  rectilinear  to  curvi¬ 
linear  motion,  then 


On  the  other  hand,  on  the  basis  of  Albrlng 's  assumption  that 
the  center  of  pressure  Is  unaltered  by  the  addition  of  fins 
or  skegs,  his  formula  for  N'^,  for  the  body  with  fins  should 
hold  for  the  bare  hull: 
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where  -7-  is  Implied  by  Albrlng  to  be  half  the  prismatic  coef- 


in  this  case  would  be  the  value  for 


TT 

flclent  of  the  body, 
the  hull  alone,  Li  ,  since  the  effects  of  the  appendages  are 

taken  as  simply  additive. 


Equations  3  are  derived  from  simplified  flow  theory 
with  the  necessary  condition  that  In  viscid  flow  the  drag  Is 
not  zero, and  so  at  an  angle  of  attack  the  lift  force  Is  not 
zero.  Under  this  condition  there  Is  an  additional  moment  due 
to  the  action  of  the  lift  force  at  a  center  of  pressure  dis¬ 
placed  from  the  CQ  by  0.  In  rotation, the  angle  of  attack 

at  the  center  of  pressure  Is  changed  by  an  amount 


1 

tan  ( — «  — g-  » - 

where 

r  =  angular  velocity 

-tr  _  -t 

""  =  IT  K 

U  =  forward  velocity 

I  =  length  of  the  hull 

R  =  radius  of  turn 

The  additional  force  coefficient  resulting  from  this  change 
In  angle  of  attack  Is 

X  r* 

and  therefore  the  rotary  force  derivative  of  eq.  1  has  an 
additional  term 


which  Is  equivalent  to  the  negative  of  the  additional  static 
moment  derivative. 
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In  treating  the  long  body  equipped  with  fins,  skegs 
or  control  surfaces,  simplified  theory  assumes  that  there  is 
no  interference  between  the  body  and  these  surfaces  so  that 
their  separate  effects  are  additive.  Tail-surface  ( skeg) 
effects  would  be  derived  from  the  lift  on  the  surface 
( obtained  either  by  measurements  or  by  assuming  the  surface 
to  be  a  wing)  by  following  the  reasoning  of  the  previous 
paragraph.  Then 

LI  assumed  or  measured 
Pt 


“  T 


(4) 


where  is  the  algebraic  distance  between  the  center  of 
gravity  and  the  center  of  pressure  on  the  tall  (x^  is 
negative) . 

Equations  3  and  4  are  the  basis  for  the  calculation 
method  used  In  the  present  report.  The  assumed  lift  rates  and 
the  other  force  and  moment  rates  derived  from  them  will  be 
discussed  in  the  following  sections. 


The  Assumed  Static  Force  And  Moment 
Rates  For  The  Bare  Hull 

The  derivative  of  the  lift  coefficient  with  respect 

to  )3,  L*  ,  Is  assumed  as  given  by  Jones'  formula  for  a  low 

% 

aspect-ratio  wing  of  span  equal  to  twice  the  ship  draft.  The 
dimensionless  lift  rate  per  unit  wing  area 
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_7r.„_7r 

Tf  ~  2 


(^) 


Is  derived  from  the  consideration  of  elliptic  load  distribu¬ 
tions  along  the  chord  and  the  span  of  a  thin  foil.  Nondlmen- 
slonallzed  on  the  basis  of  the  area  H  x  -t,  the  static  rates 
of  eq.  3  become 


ttH 

T 


ttH 


% 


>  (5) 


where  is  the  drag  coefficient  obtained  from  the  Taylor 
Standard  Series  curves  of  resistance and  x^  is  taken  as  the 
distance  of  the  center  of  area  of  the  bare  hull  profile  from 
the  center  of  gravity. 

Tsakonas'*  and  Martin^  found  the  Jones  formula  to  be 
a  good  approximation  of  the  static  lift  rate  of  the  Standard 
Series  hulls.  The  fact  that  the  formula  does  not  predict  the 
values  obtained  from  measurements  on  flat  plates  of  the  same 
profile  does  not  detract  from  its  usefulness  for  cambered 
wings.  Blspllnghoff  proves  that  an  elllptlcally  loaded  wing 
must  have  an  elliptic  planform  and  he  points  out  that  it  is 
the  so-called  flat-plate  chordwise  distribution  which  "compares 
very  favorably  with  the  measured  chordwise  distribution  of 
pressure  difference  over  a  slightly  inclined,  thin,  uncambered 
airfoil  developing  the  same  lift.  The  only  significant  dis¬ 
crepancies  come  from  within  a  few  per  cent  chord  lengths  of 
the  singularity  [at  the  leading  edge]."^^  Tsakonas'*  referred 
to  Crabtree’s^®  work  in  substantiation  of  this  statement. 
Crabtree  had  found  that  the  pressure  distribution  over  a  thin 
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plate  ( less  than  12^  thickness  to  chord  ratio)  at  various 
incidences  showed  a  pronounced  suction  peak  near  the  leading 
edge  with  a  consequent  steep  adverse  pressure  gradient  and 
laminar  boundary  layer  separation.  Since  the  size  and  form 
of  the  separation  region  or  "bubble"  has  a  large  effect  on 
the  lift,  Tsakonas  cautioned  against  use  of  the  experimental 
measurements  on  flat  plates  to  predict  forces  in  the  case  of 
the  ship- wing  analogy. 


The  Assumed  Rotary  Force  And  Moment 
Rates  For  The  Bare  Hull 

From  simplified  flow  theory  the  rotary  force  de¬ 
rivative  for  a  body  of  revolution  without  appendages  is  zero 
if  the  body  has  fore-and-aft  symmetry.  For  a  bare  ship  hull 
without  fore-and-aft  symmetry,  on  the  other  hand,  eq.  3  shows 
that; 


Ki  ■  -  if-  •  "h 


- 

1 


H  ^ 

and  the  measured  force  derivative  with  respect  to  r'  (eq.  2)' 


7^  =  -  (m*  +  m»)  - 


L' 


where  Xp  is  different  from  zero  although  small. 
From  eq.  3a J 


which  is  very  small.  In  addition,  as  shown  by  Martin,^  there 
is  a  term  due  to  the  asymmetry  of  the  rotary  added  mass  of  en¬ 
trained  water,  m*.  The  moment  due  to  this  added  mass  is 
z 

-m^r’  x./l  where  x  is  the  distance  from  the  center  of  gravity 
of  the  hull  to  the  CG  of  the  lateral  added  mass. 
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(6) 


Then  eq.  3a  becomes 


-  -•‘iz-  \  (6*) 

If  eq.  3b  Is  used 

"r-  -  -KI-  \ 

where  x.^1  Is  half  the  prismatic  coefficient.  The  values  ob¬ 
tained  from  eq.  6b  have  been  found  to  be  much  closer  to  ex¬ 
perimental  results  than  those  calculated  by  eq.  6a.  The 
second  term  of  eq.  6b  is  approximately  I-I/2  times  the  first. 
For  this  reason,  eq.  6b  has  been  taken  as  the  assumed  rotary 

moment  rate.  As  in  ref.  1,  m*  and  x  have  been  estimated  in 

z 

the  following  way: 


m^  =  kg  ^ 


C  h^  dx 
s 


kg^k*  are  Lamb’s  coefficients  of  accession  to  Inertia, 
lateral  and  rotational. 

Xg,Xb  are  the  x  coordinates  of  the  stern,  bow. 
h  =  local  draft 

C  =  two-dimensional  lateral  added  mass  coefficient. 

It  is  determined  at  each  section  from  the  curves 
on  two-dimensional  forms  of  Lewis’  sections  by 
Prohaska. 
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X  =  LOG 


C„h^xdx 

s 


C„h^dx 

O 


The  Stability  Derivatives 
For  Hulls  With  Skegs 


The  force  and  moment  derivatives  for  the  hulls  with 
skegs  are  taken  as  a  simple  addition  of  the  corresponding 
forms  of  eqs.  6  and  4: 

Li  =  Li  +  Li  \ 

^  %  ^T 


Yi  =  Y'  +  L) 


X 


T? 


=  Yi, 


^rp 

~r  \ 


.^T  a 


(7) 


n;.  =  n;.  -  (i^)  J 

In  the  present  nondimenslonallzed  form  of  Jones*  formula 


h^  =  maximum  skeg  height 

The  center  of  pressure  of  the  skeg  will  be  assumed  at  the 
after  end  of  the  skeg.  Since  the  length  of  the  skeg  is  small 
In  comparison  with  the  length  of  the  hull,  any  error  involved 
in  this  assumption  will  be  small. 
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For  the  case  of  models  with  skegs  extending  to  the 
stern  (Pig.  2),  these  formulas  give  essentially  the  same  re¬ 
sults  as  those  obtained  by  the  method  suggested  by  Martin.^ 

He  Included  the  following  terms  In  his  force  and  moment 
equations,  respectively,  to  account  for  the  skeg  effect: 

-  m(x„)u^  +  k’r  x„)  (from  eq.  2  of  ref.  1  for  P  ) 

s  X  «  y  3  y 

(9a) 

and 

-  X  m(x  )u  (k  u  +  k’r  x  )  (from  eq.  3  of  ref.  1  for  N) 

s  s  X  y  s 

(9b) 

where 

ni(x  )=  two-dimensional  lateral  added  mass  at  the  stern 
®  per  ft 

u  =  U  cos  ^  ^  V 
u  -U  Gin  iS  ^  -U^ 

•X 

r  =  angular  velocity  of  the  ship. 

Vlien  the  section  at  the  stern  Is  that  of  a  flat  plate, 

(see  Kuerti,  McFadden  and  Shanks^®) 

and  this  term  is  numerically  the  same  as  that  given  hy  the 
Jones  formula  for  a  low  aspect-ratio  wing  of  span  equal  to 
twice  the  ship  draft.  After  nondimensionalizing,  since 
^  1^  eq.  9a  for  the  added  force  is  approximately 

\  t) 

and  eq.  9b  for  the  added  moment  Is 
/  ^T\ 

-f) 
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where  Li  Is  given  by  eq.  8.  Then  the  force  and  moment  rates 

pt 

are  those  given  in  eq.  4. 

The  Stability  Indices 

The  criteria  for  Inherent  dynamic  stability  of  a 
free  body  moving  on  straight  course  in  the  horizontal  plane 
are  the  damping  exponents  and  Og  in  the  solution 


OiS 


02  3 


0lS 


^  r'  =  r^  e 


+  r^  e 


02  S 


jS  =  e  +  i32  e 
of  the  homogeneous  linearized  equations  of  motion^® 


( 


m* 


Y;,)  r»  -  m^  -  Y'  ^  =  0 


y  "^s 


"z  -  N^.r'  -  |3  =  0 


(10) 


where  Y^,  and  N^,  are  defined  in  eqs.  6  and  6b  for 

the  bare  hull  and  in  eq.  7  for  the  hull  with  skegs, 

=  m^  +  m{ 


"ly  =  ">6  +  '"a 


n'  = 


0  %  , 

■jg*  (assuming  the  radius  of  gyration  is  equal  to 


H 


Iz  =  k- 


^1 

¥  / 


Cg  h^x^dx  (ref.  l) 


and 


^.c 


s  =  5i“' 


9r  ’  Ut 

,,  s  =  -j- 


s  3s  J 

Solution  of  the  characteristic  equation  of  eq.  10  gives  the  roots 
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'l^2" 


■(  n«Y « -m^N « , ) ±  ^  ( n » Y « -m^N » . )  ^+4n «m «  [N^ . Y^+(  m^-Y^ , ) N »  ] 


2n' 

z 


m; 


(11) 


If  (Ti  and  (Ja  ( or  their  real  parts)  are  both  negative, 
the  motion  Is  Inherently  stable  In  that  an  Initial  disturbance 
damps  out  exponentially;  the  more  negative  the  exponents,  the 
sooner  It  damps  out  and  hence  the  greater  the  stability.  If 
Oi  Is  positive,  the  motion  Is  unstable  and  the  hull  cannot 
keep  to  a  straight  course  without  application  of  a  corrective 
rudder . 


The  turning  characteristics  of  a  hull  In  turns  that 
are  not  too  tight  (when  nonllnearltles  can  be  neglected)  can 
be  predicted  qualitatively  from  the  Index.  A  more  dynam¬ 
ically  stable  hull  will  turn  In  a  larger  radius  under  a  given 
applied  rudder  force  than  will  a  less  stable  hull.  Conversely, 
the  more  stable  hull  will  require  greater  rudder  force  than 
the  less  stable  hull  to  turn  In  a  given  radius.  On  the  other 
hand,  an  unstable  body  may  turn  In  the  opposite  direction  to 
that  Indicated  by  the  applied  rudder  and  will  need  a  com¬ 
paratively  large  force  to  bring  It  around  (until  It  becomes 
stable  In  the  turn,  but  stability  In  turn  Is  outside  the  scope 
of  this  report) . 
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PRESENTATION  AND  DISCUSSION  OF  RESULTS 


The  experimental  data  measured  In  1951 

(unpublished),  and  1959^'*  are  shovm  In  Pigs.  A-1  to  A- 48  In 
the  Appendix,  plotted  In  dimensionless  form  versus  dimension¬ 
less  angular  velocity  r*  =  -t/R  and  yaw  angle  jS.  Figures 
A-1  to  A-19  present  the  measured  lateral  force  coefficients 

P’  =  -(m'  +  m|)  r’  +  y*  and  Figs.  A-20  to  A-38  the  measured 
y  o 

yawing  moment  coefficients  N’  for  the  bare  hulls.  Figs.  A-39 
to  A-48  show  these  quantities  for  three  of  the  hulls  equipped 
with  skegs  but  no  other  appendages.  Skeg  20  or  A  Is  the  full 
skeg  extending  to  the  after  perpendicular  (Sta.  20),  skeg  18 
or  B  extends  to  Sta.  18,  0.1  length  from  the  after  perpendi¬ 
cular,  and  skeg  I7  or  C  to  Sta.  17,  0.15  length  from  the  after 
perpendicular. 

The  experiments  had  been  conducted  at  speeds  ranging 
from  length-Proude  numbers  of  0.I6  to  0.23,  for  which  range 
the  Implication  of  the  analytical  method  that  wave  making  may 
be  neglected  is  valid.  In  that  speed  range  the  ratio  of  force 
or  moment  to  the  square  of  speed  changes  only  slightly  so  that 
data  measured  at  different  speeds  within  the  range  can  be  com¬ 
pared  legitimately. 

Table  I  gives  the  pertinent  characteristics  of  the 
eight  hulls  and  the  necessary  Information  for  the  calculations. 
The  resistance  coefficients  D^  =  2  (r^  4-  R^)/pU^tH  were  taken 
from  the  Taylor  Standard  Series  curves^'*  as  the  average  of  the 
slightly  different  values  for  Proude  numbers  O.16  and  0.23 . 

The  calculated  Y^/S  and  N^jS,  force  and  moment  coef¬ 
ficients  respectively  at  r*  =0,  are  also  shown  on  the  data 
charts  in  the  Appendix.  The  data  plotted  at  various  /3  versus 
r’  are  faired  to  the  calculated  values  at  r'  =0  with  no 
stretch  of  the  Imagination.  The  static  force  and  moment  rate 
coefficients  are  the  same  whether  predicted  by  the  einalytlcal 
method  or  by  rotating-arm  data. 
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Figures  4  through  9  In  the  text  are  summary  charts 
comparing  experimentally  obtained  values  for  the  rotary  as 
well  as  static  derivatives  with  those  calculated  by  the 
analytical  method.  These  values  are  shown  for  varying  length- 
draft  -t/H  and  beaim-draft  b/h  ratios. 

The  calculated  SP^/Sr  and  stability  Indices 

0i,2  predict  the  variations  with  -t/H  and  B/H  correctly  but 
underestimate  the  actual  magnitudes  slightly  In  the  bare-hull 
cases,  more  In  the  cases  of  hulls  with  skegs.  The  quantita¬ 
tive  predictions,  however,  are  on  the  conservative  side. 

These  comparisons  suggest  that  the  simple  method  adopted  here 
can  be  useful  for  estimating  the  stability  of  a  given  vessel. 
The  good  results  Justify  use  of  the  shlp-wlng  analogy  and  the 
Jones  formula  for  the  lift  on  the  ship  as  a  low  aspect-ratio 
wing. 

The  following  are  specific  deductions  from  the 

charts: 


1.  The  static  force  derivative  varies  Inversely 
with  -t/H  (or  directly  with  aspect  ratio)  and  directly  with 
B/h  although  B/h  has  very  slight  effect.  Is  Increased  by 
adding  skeg  area.  An  Increase  In  Y^  Is  In  the  direction  of 
greater  stability. 

2.  The  static  moment  derivative  also  varies  in¬ 
versely  with  -t/H  eind  directly  with  b/h.  This  destabllzlng 
moment  rate  Is  reduced  mainly  by  Increasing  -t/H  and  by  adding 
skeg  area  at  the  stern. 

3.  The  rotary  force  derivative  dF^/dr’  becomes 
less  negative  or  more  positive  by  Increasing  -t/H,  by  de¬ 
creasing  b/h  and  by  Increasing  skeg  area  at  the  stern.  As  In 
the  case  of  N^,  Increasing  -t/H  or  decreasing  B/h  Is  In  the 
direction  of  greater  stability.  The  variation  with  B/h  arises 
from  the  variation  In  the  longitudinal  virtual  mass  coeffl- 


R-945 

-16- 


dent.  This  is  the  most  Important  effect  for  the  bare  hulls, 
and  is  as  important  as  the  skeg  effect  in  the  case  of  hull 
with  full  skeg. 

4.  The  rotary  moment  derivative  N\  is  independ¬ 
ent  of  B/h.  It  becomes  less  negative  with  increasing  t/H 
and  more  negative  ( towards  greater  stability)  when  skeg  area 
is  added  at  the  stern.  The  variation  with  t/H  is  slight  in 
the  case  of  the  bare  hulls,  more  pronounced  for  the  hulls 
with  skegs. 

5.  The  stability  indices  which  combine  the 

effects  of  static  and  rotary  force  and  moment  rates,  show 
that  stability  depends  almost  entirely  on  B/h  and  very  little 
on  i/h  in  the  range  tested.  Stability  Increases  as  B/h  de¬ 
creases.  On  the  other  hand,  an  Increase  in  b/h  would  result 
in  greater  turning  ability. 

Figure  10  is  a  summary  chart  comparing  the  calcu¬ 
lated  Y^,  N^,  Y^i,  and  Oi  with  the  values  from  measure¬ 
ments  obtained  by  Tsakonas^  in  1959  on  the  842  hull,  without 
skeg  and  with  three  skegs  of  different  sizes.  The  calculated 
static  derivatives  and  are  Identical  with  those  obtained 
from  experimental  results.  This  Is  also  shown  in  Figs.  4  and 
7.  Figure  10  shows  that  while  the  calculated  and  experimental 
magnitudes  of  the  rotary  derivatives  and  stability  index  differ 
slightly,  the  stability  predictions  are  conservative  estimates. 
It  is  also  seen  that  the  analytical  method  can  predict  the 
trend  in  stability  with  Increase  in  skeg  area.  Figures  11  and 
12  for  the  hitherto  unpublished  1951  data  confirm  these  con¬ 
clusions  . 
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CONCLUSIONS  AND  RECOMMENDATIONS 


An  analytical  method  for  estimating  course  stability 
and  turning  qualities  of  ships  has  been  developed  and  com¬ 
pared  with  available  experimental  data  on  a  series  of  eight 
hulls,  the  840  Series,  of  the  same  length  and  prismatic  but 
varying  In  draft  and/or  beam  and  displacement.  The  analyti¬ 
cal  method  combines  simplified  flow  theory  with  low  aspect- 
ratio  wing  theory  and  makes  use  of  Albrlng’s  empirical  modi¬ 
fication  for  the  rotary  moment  derivative. 

Encouragingly  good  correlation  Is  shown  between  the 
calculations  by  this  method  and  the  results  based  on  the  ex¬ 
perimental  data.  However,  since  the  840  Series  Is  a  family 
of  hulls  of  the  same  prismatic,  Albrlng’s  modification  for 
the  rotary  moment  rate  which  is  a  function  of  prismatic  coef¬ 
ficient  has  not  been  fully  tested.  Necessary  refinements  of 
the  method  must  wait  on  analysis  of  data  on  hulls  of  other 
prismatic,  with  and  without  skegs  or  deadwood  aft. 

It  Is  recommended  therefore  that  all  available  data 
on  hydrodyneimlc  forces  and  moments.  In  turn,  for  other  hull 
forms  with  different  prismatic  as  well  as  length-draft  and 
beam-draft  ratios  be  assembled  and  analyzed  with  a  view  to 
checking  and  refining  the  method.  It  would  also  be  advisable 
to  do  further  experimentation  In  still  water  on  families  of 
hulls  such  as  the  Series  60.  The  latter  series  has  been  the 
subject  of  extensive  tests  to  determine  resistance,  bending 
moments  and  sea-keeping  qualities  In  waves  but  the  stability 
on  straight  course  and  behavior  in  turn  have  so  far  been  in¬ 
vestigated  for  only  one  of  the  series.^® 
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FIGURE  8.  840  SERIES  HULLS  WITH  FULL  SKEG  (20) 
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FIGURE  1 1 .  COMPARISON  OF  CALCULATED  AND  EXPERIMENTAL 

STABILITY  DERIVATIVES  AND  INDICES  FOR  646  HULL  WITH 
OUT  AND  WITH  SKE6  (1951  DATA) 
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